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Abstract. We here discuss the process of opinion formation in an open community where agents are made 
to interact and consequently update their beliefs. New actors (birth) are assumed to replace individuals 
that abandon the community (deaths). This dynamics is simulated in the framework of a simplified model 
that accounts for mutual affinity between agents. A rich phenomenology is presented and discussed with 
reference to the original (closed group) setting. Numerical findings are supported by analytical calculations. 

PACS. 87.23. Ge Dynamics of social systems - 05.45.-a Nonlinear dynamics and nonlinear dynamical 
systems 

1 Introduction cial groups. Several models have been proposed in the past 

to reproduce the key elements that supposedly drive the 

„ . . , . , ,. , , ,1 • n 1 1 process of opinion making [Tl[21 . The problem of providina' 

Upmion dynamics modeling represents a challenging held r i o s— i.— i 

, ., r . ,• ,• 1 , • 1 T ■ an adequate experimental backup to such developments is 

where ideas trom statistical physics and non-linear science 

, ., , T 1 i 1 i 1 ii r indeed an open one. and more work is certainly needed to 

can be possibly applied to understand the emergence or 

,, ,. 1 , . ,., ..... eventually assess the interpretative ability of the proposed 

collective behaviors, like consensus or polarization m so- ■^ ^ j t- t- 
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mathematical formulations, following e.g. the guidelines of To this end, we refer to the model presented in [7], 

[3ll4] . Models based on interacting agents display however where the role of affinity among individuals is introduced 

a rich and intriguing dynamics which deserves to be fully as an additional ingredient. This novel quantity measures 

unraveled. the degree of inter-personal intimacy and sharing, an ef- 



fect of paramount importance in real social system [8] . 
Indeed, the outcome of an hypothetic binary interaction 
relies on the difference of opinions, previously postulated, 
but also on the quality of the mutual relationships. The 
affinity is dynamically coupled to the opinion, and, in this 
respect, it introduces a memory bias into the system: affin- 
ity between agents increases when their opinion tends to 
converge. 

The aforementioned model is here modified to accom- 



Opinion dynamics models can be classified in two large 
groups. On the one hand, opinions are represented as dis- 
crete (spin-like) variables where the system behaves simi- 
larly to spin glasses models [5j. On the other, each individ- 
ual bears a continuous opinion which span a pre-assigned 
range [5]. More recently, a new framework for a discrete 
but unbounded number of opinions was proposed in [2] 

and shown to nicely complement the picture. In all the 

, 1 1 1 i • 11 1 modate for a death/birth like process. In this formulation, 

above approaches, a closed system is generally assumed, ' 

^1 ^ ^1 1 r i ■ 1 i ■ i ^ M agents are randomly eliminated from the system, ev- 

meanmg that the same pool of actors is made to interact ^ j ^ 

, . ^1 , , . mi • 1 • i ill cry T time steps. When an agent exits from the commu- 

durmg the evolution, ihis can be interpreted by assum- j r o 

,1 , ii . 1 i 1 • 1 ,■ 1 / ■ • nity (virtually, dies), he is immediately replaced by a new 

mg that the inherent dynamical timescales (e.g. opinion ■' ^ ji /; j t- j 

, . N 1 r X ii ii ■ i 1 element, whose opinion and affinity with respect to the 

convergence time) are much faster than those associated ^ t- j ±- 

, ,1 / ■ i- 1 ■ ii / 1 j_i \ -11 group are randomly assigned. As we shall see, the pertur- 

to the processes (e.g. migration, birth/death) responsible <= ^ j t, i t- 

r I.J, , . r ii ... ,1 1 ^^ bation here prescribed alters dramatically the behavior of 

tor a modification ot the group composition, these latter 

rr ^ T ■ 1- ii r r liioi the system, with reference to the ideal close-system config- 

ettects having being therefore so far neglected. Such an -^ jo 

T ., i- ■ i • 1 i 1 ii J u i uration. To understand such modifications via combined 

implicit assumption is certainly correct when the debate 

1 1 , 11 -i r ■ 1- • 1 1 ii 1 numerical and analytical tools, constitutes the object of 

bound to a small community ot individuals, thus mak- 



is 



the investigations here reported. 



ing it possible to eventually achieve a rapid convergence 

towards the final configuration. Conversely, it might prove We shall mainly explore the parameters setting that 

inaccurate when applied to a large ensemble of interacting would lead to an asymptotic consensus state (all agents 

agents, as the process becomes considerably slower and ex- eventually bearing the shared opinion 0.5), in absence of 

ternal perturbations need to be accounted for. Given the the external perturbation which is here object of investiga- 

above, it is therefore of interest to elucidate the open sys- tions. The underlying model [7j displays however a richer 

tem setting, where the population is periodically renewed, phenomenology, exhibiting in particular stable polarized 
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states in the late time evolution and long-lived metastable superposition of the so-called social distance, dij, and a 

regimes. Though it would be extremely interesting to ex- stochastic contribution r]j, namely: 

tend the present analysis and hence cover those additional „„ ,* , ^ 

Dlj^4j+ViiO,c7). (1) 

scenarios, we chose to only briefly touch upon this issue 
when commenting the details of the transition between 
single and fragmented phases. 

The paper is structured as follows. We first introduce 
the model, then we present the obtained analytical and 
numerical results and, finally, we sum up and draw our "ij = ^Oyd — a^j) j = 1, ...,N j t" i , (2) 



Here ?/j(0,cr) represents a normally distributed noise, 
with mean zero and variance a, the latter being named 
social temperature. The social distance is instead defined 
as: 



conclusions. 



2 The model 



withZ\a*j =0*-0*. 

The smaller the value of d*- the closer the agent j 
to i, both in term of affinity and opinion. The additive 



T ;i r 1, . .1, • ;i 11 • 1 • noise 77,(0, cr) acts therefore on a fictitious ID manifold, 

m the toilowmg, we will review the model previously m- -^ ^ ' 

, , . 1=, , , ,, , ,.,. 1 r , ,1 , which is introduced to define the pseudo-particle (agent) 

troduced m [7j and present the additional features that \ o / 

, , . . • mi • ; ; 1 1 .1 interaction on the basis of a nearest neighbors selection 

are here under inspection, i he interested reader can thus 

r , ,, • • 1 pti r jj-i- 1 i mechanism and, in this respect, set the degree of mixing 

refer to the original paper [7j tor an additional account on ^ i- t o o 

,1111,.,. in the community, 

the model characteristics. 

When the two agents i and j are extracted on the basis 



Consider a population of N agents and assume that at 
time i they bear a scalar opinion O- e [0, f]. We also intro- * * . ~ . _ 



of the recipe prescribed above, they interact and update 

1 .lATAT,- 1 1, ;■ t 1 1 , their characteristics according to the following scheme I 

duce the N xN time dependent matrix a , whose elements 



(3) 



alj ' = 0L% + a% [1 - ay 1^2 {AO^) , 



where the functions A and i^2 respectively read: 

A (a*,) = e {a% - ae) (4) 



a*- belong to the interval [0, f]. The quantities a* specify O*"*"^ = 0\ — \A0\ Fx (a* ) 

the affinity of individual i vs. j, at time t: Larger values | ^,t+i _ ^,t , ^,t r-, ^,t 

of a*j are associated to more trustable relationships. 

Both the affinity matrix and the agents opinions are 
randomly initialized time i = 0. At each time step t, two 
agents, say i and j, are chosen according to the following ^'^"^ 

extraction rule: first the agent i is randomly identified, j^ {AO ■'] = O (AO — \AO^\) — i\AO^-\ — AO ) (5) 
with a uniform probability. Then, the agent j which is ^ The evolution of the quantities Oj(i) and Qji(t) is straight- 
closer to i in term of the social metric D^- is selected forwardly obtained by switching the labels i and j in the equa- 
for interaction. The quantity D^' results from the linear tions. 
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and the symbol 0{-) stands for the Heaviside step-functionlj 
More specifically, A is or 1 while /^2 is —1 or 1, de- 
pending on the value of their respective arguments. In the 
above expressions ac and AOc are constant parameters. 
Notice that, for etc — > 0, the opinion is formally decoupled 
from affinity in Q being Fi = I irrespectively of the ac- 
tual value of a*, , and the former evolves following the Def- 
fuant et al. scheme [6] with convergence rate fj. = 0.5 and 
interaction threshold d = I (confidence bound). The lat- 
ter scheme pioneered the broad class of models inspired to 
the so-called bounded confidence hypothesis, an assump- 
tion which, though revisited, also enters the self-consistent 
scenario of Eqs.Q. 

In [7] , a preliminary analysis of the qualitative behav- 
ior of the model as a function of the involved parameters is 
reported. Asymptotic clusters of opinion are formed, each 
agglomeration being different in size and centered around 
distinct opinion values. Individuals sharing the same be- 
lieves are also characterized by a large affinity scores, as 
it is exemplified in Fig. [TJ 

More quantitatively, the system is shown to undergo 
a continuous phase transition: above a critical value of 
the control parameter {aac}^^''^ it fragments into several 
opinion clusters, otherwise convergence to a single group is 




°04 




Fig. 1. Upper panel: Typical evolution of the opinion versus 
time, i.e. number of iterations and the asymptotic distribution 
of agents' opinion in a closed community. Bottom panel: Final 
affinity matrix, represented by the underlying network: node i 
is linked to node j if aij > a = 0.8. Here a — A- 10^'', AOc = 
0.5, Uc = 0.5 and p — 0, i.e. no agent can leave the group. 
Initial opinion are (random) uniformly distributed within the 
interval [0,1], while a''ij is initialized with uniform (random) 
values between an 0.5. 



.a 



numerically shown to occur |j. We shall here simply notice 



^ In [7], the switchers A and /2 are smooth functions con- 
structed from the hyperbolic tangent. We shall here limit the 
discussion to considering the Heaviside approximation, which 
is recovered by formally sending /3i,2 to infinity in Eqs. (3) and 

(4) of m- 



■^ Strictly speaking, it should be noted that the fragmented 
state is metastable, if the mean separation between the adja- 
cent peaks is smaller than the interaction distance AOc- There 
always exists in fact a finite, though small, probability of se- 
lecting two individuals which belong to different agglomera- 
tions. When the above condition applies, i.e. when the agents' 
opinions are closer than the threshold amount AOc, such rare 
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that a significant degree of mixing (large social tempera- 3 Results 

ture a) brings the system towards the single-cluster final 

n , • Numerical simulations are performed for a system of A'^ = 

100 individuals and its evolution monitored Q Qualita- 
Starting from this setting, we introduce the birth/death ^i^^ly, the system shows the typical critical behavior as 
process, which in turn amounts to place the system in con- observed in the original formulation 0. However, pecu- 
tact with an external reservoir. The perturbation here hy- ^'^"^ distinctions are found, some of those being addressed 
pothesized is periodic and leaves the total number of agent '^ ^he following discussion. First, an apparently smooth 
unchanged. Every T time steps (i.e. encounter events) M transition is also observed within this novel formulation, 
agents, randomly selected, are forced to abandon the sys- ^^ich divides the mono- and multi-clustered phases. In- 
tem (death). Every removed individual is instantaneously terestingly, the transition point is now sensitive to the de- 
replaced by a new element, whose initial opinion and affin- P^^^ture density p. To further elucidate this point, we draw 
ity are randomly fished, with uniform probability, from the ^^ ^ig. [2] the average number of observed clusters versus 
respective intervals [0, 1] and [0, a^.^J- Further, we intro- ^ rescaled temperature. A clear transition towards an or- 
duce p = f to characterize the departure density, a crucial ^cred (single-clustered) phase is observed, as the temper- 
quantity that will play the role of the control parameter in ^ture increases. The parameter a, in Fig. m plays the role 
our subsequent developments. As a final remark, it should "^ ^^ effective temperature, and it is numerically adjusted 
be emphasized that no aging mechanisms are introduced: ^o make distinct curves collapse onto the same profile, 
agents are mature enough to experience peer to peer en- ^^ich hence applies to all values of p. The inset of Fig. M 
counters from the time they enter the system. ^hows that there is a linear correlation between a, and p. 



* The chosen value of N could be in principle considered 
too small to allow us extracting sound statistical information 
from the model at hand. As we shall however discuss, already 
at such relatively small value of A'^, one observes a satisfy- 



encounters produce a gradual increase of the mutual affinity ing matching between numerics and statistical based predic- 
scores, a tendency which asymptotically drives a merging of tions. No substantial differences are detected when simulating 
the segregated clusters, as ruled by Eqs. Q. This final con- a larger system, this observation motivating our choice to stick 
vergence is eventually achieved on extremely long time scales, to the N = 100 case study. Notice also that potentially inter- 
diverging with the number of agents. Socially relevant dynam- esting applications in social sciences would often deal with a 
ics are hence bound to the metastable phases, which are being finite, possibly limited, number of agents, as for the case being 
investigated in Fig. [2] addressed at present. 



O o = 0.29 p + a 



'1? 
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The larger the departure density p, the larger the effcc- , — . — . — . — . 

, xlO^ 

10 - ,, 7r 

6 

tive temperature CTc- In other words, when p is made to ■ " +, 

increase (i.e. the system is experiencing the effect of a more 
pronounced external perturbation), one needs to augment ^ 

the degree of mixing, here controlled by the social tem- 
perature a, if a convergence to the final mono-cluster is 

10°l ' ^ 

sought. The death/birth process is in fact acting against 

the thermal contribution, which brings into contact oth- 

Fig. 2. Main plot: Average number of clusters as function 
erwise socially distant individuals. While this latter effect 

of the rescaled quantity -2- for different values of the density 
enhances the chances of convergence, the former favors the 

p. Simulations have been performed with parameter values: 
opposite tendency to spread. 

p = (O), p = 0.0025 (*), p = 0.0033 (+), p = 0.005 (x). 

To further elucidate the role of the external pertur- ^ ^ q g^ („^ ^^^ ^ ^ q 05 (^^ j^ ^H simulations, here and after 

bation, we shall refer to the dynamical regime where the reported, unless otherwise specified, a, = AO, = 0.5; O^ and 

agents converge to a single cluster. When p is set to zero, a°j are random variables uniformly distributed in the intervals 

the final shared opinion is 0.5 to which aU agents eventu- [0, f] and [0,amaa:l - being a„aa: = 0.5 - respectively. Inset : 

ally agree, see Fig. [3^. In other words, the final distribu- Cc as a function of p. The open circles, o, represent the values 

tion is a Dirac delta, with the peak positioned in O = 0.5. calculated from the transitions shown in the main plot. The 

Conversely, for positive, but smaU, values of p, the final dotted lines represent the best linear fit : a^ = 0.3p + aco, with 

distribution of opinions presents a clear spreading around ^co = 5.5 • 10"* . Notice that a, = a,„ is eventually recovered in 

the most probable value, stiU found to be 0.5. This see- ^^^^ closed-system setting, which in turn corresponds to p = 0. 
nario is clearly depicted in Fig. [SJd. For larger p, when 

the birth-death perturbation becomes more frequent, the t^ hi r ^i ^ i _l _li ^ 

_tor small values oi the control parameter p, the stan- 

opinion profile cannot relax away from the initial distri- t , , ■ ^- r ii i ^ i i- n ^ 

>: >^ -^ dard deviation v of the cluster scales proportionally to 

bution, the agents believes being uniformly scattered over ^ „. ^^ ^- • • i- ^ ^i ^ ^i ^- ^■l_ 

" o J ^^ ^ jg |3j jxiunierics indicate that the proportionality co- 

the allowed interval, i.e. [0,11. xc • ^ j. n t ^i -^ t -i 

<- ^ ' eihcient gets smaller, as ac grows, in the opposite hmit. 

The associated standard deviation v is deduced, from a namely for large values of the density p, the standard de- 
series of simulations, and shown to depend on the selected viation v rapidly saturates to a asymptotic value, Vc- The 
value of p. The result of the analysis is reported in Fig. 31 latter is universal, meaning that it neither scales with ac, 
where the calculated value of u (symbols) is plotted versus nor it does with the social temperature a. Our best nu- 
the departure density amount p. merical estimates returns, Vc — 0.28 ~ l/\/T2 which, as 
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Fig. 3. Opinion evolution versus time, the latter being quantified through the number of iterations. The black stars represent 
the mean opinion and the error bar the standard deviation of the opinion distribution. In the insets, the histogram of asymptotic 
distribution of agents' opinion. In panel a) the system is closed (i.e. p = 0), in b) p = 5 • 10~^. 



expected, corresponds to the standard deviation of the 
uniform distribution in the interval of length 1. 
The solid lines in Fig. |4] represent the function : 

,2 . ^ (6) 

which straightforwardly follows from an analytical argu- 
ment, developed hereafter. In the above expression Tc, 
stands for an effective estimate for time of convergence 
of the opinion cluster, and is deduced via numerical fit 
(see caption of Fig. |4]and [10] for further details). In [9], 
working within the Deffuant's scheme [5], i.e. closed com- 
munity case without affinity, the time needed to form a 
coherent assembly from a sequence of binary encounters 
was shown to diverge with the population size iV, with a 
super-linear scaling. Moreover, it was also proven that the 
affinity slows down the convergence rate, a fact that can 
be successfully captured by accounting for an additional 



dependence of T^ over oic'- The larger a^ the longer the 
convergence time, as reported in [J. A comprehensive dis- 
cussion on the analytical derivation of Tc{ac) falls outside 
the scope of the present discussion and will be presented 
in a forthcoming contribution [llj . 

Before turning to discuss the analytical derivation of 
Eq. ([5]), we wish to test its predictive adequacy with ref- 
erence to the two limiting cases outlined above. Indeed, 
for p << 1 and T « Tc, Eq. ([6]) can be cast in the 
approximated form: 



Tcjac) 

24iV 



Ity/p, 



(7) 



which presents the same dependence of v versus y/p, as 
observed in the numerical experiments. Moreover, the co- 
efficient "ft is expected to decay when increasing the cutoff 
in affinity ac, in agreement with the numerics. For p >> 1, 
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Eq. ^ implies: 




1 
12 ' 



(8) 



Fig. 4. Main: The standard deviation of the final mono-cluster 
w as a function of the departure density p. Each point results 
from averaging out 100 independent runs, with A'^ = 100. Sym- 
bols refer to numerical simulation performed with different val- 
ues of Oc, M and a : (*) a^ = 0.5, M = l,o- = 0.07, (O) 
Oc = 0.5, M = 2, cr = 0.07, (□) a^ = 0.5, M = 6, a = 0.07, 
(•) Qc = 0.5, M = 2, cr = 0.25, (v) ac = 0.5, Af = 2, a = 1, 
(o) Qc = 0.3, Af = 2, a = 0.28, (A) Oc = 0.3, M = 6, cr = 0.28 
and (x) Oc — 0, M = 2, a = 0.28). The solids line refers to 
the theoretical prediction (|6]): The free parameter Tc is numer- 
ically fitted and results in Tc = 8030 for Oc = 0.5, Tc = 1886 
for Qfc — 0.3 and Tc = 1470 for Oc = 0. This values are in good 
agreement with the simulation results [lO], which confirms the 
validity of the proposed analytical scheme. Inset: v versus y/p, 
for p G [0, 0.5]. This zoomed view confirms the correctness of 
the scaling dependence derived in ((7]). 



thus returning the correct result. 

To derive Eq. ^ let us suppose that at time t the 
death/birth process takes place and the system experience 
an injection of new individuals. Label with vt the stan- 
dard deviation of the agents opinion distribution f^{0), 
at time t. It is reasonable to assume that /*(0) is cen- 
tered around 1/2. After T interactions between agents, 
when the next perturbation will occur (M agents are ran- 
domly removed from the community and replaced by Af 
new actors with random opinion and affinity scores) the 
distribution has been already modified, because of the un- 
derlying dynamical mechanism specified through Eqs. ([3]). 
More concretely, the opinions slightly converge around the 
peak value 1/2, an effect that certainly translates into a re- 
duction of the associated standard deviation. To provide a 
quantitative estimate of such phenomenon, we recall that 
in the relevant (O, t) plan, the convergence process fills an 
ideal triangular pattern, whose height measures Tc- This 
topological observation enables us to put forward the fol- 
lowing linear ansatz: 



^t+T 



Vt 1 



T 

t: 



(9) 



where uf?^" labels the standard deviation just before the 

. . n 

insertion of the next pool of incoming agents |2| 



^ Numerical simulations (not reported here) show that in the 
closed model, the standard deviation of the opinions' distribu- 
tion exhibits a exponential decay as a function of a power of 
time. This latter is approximately interpolated by the proposed 
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Recalling that the newly inserted elements arc uni- that the final fate of the system is not affected by the 
formly distributed, and labeling with G(-, •) the opinion time when the perturbation is first applied, t.i„. This con- 
distribution (the two entries referring respectively to mean elusion is also confirmed by direct numerical inspection: 
and the standard deviation), the updated variance is: The asymptotic value of the standard deviation v does 



M f^ f ^ l\^ „ J^ot depend on ti„, but solely on p. Even in the extreme 

condition, when the death/birth perturbation is switched 



"--=ivX 1^-2' '^ 



N - M /"^ „ /I _„,A /„ 1^^ 



/^ _ conv \ I (^ _ _ \ il(~\ 

N Jq \ 2 ' *^"^ J \ 2 J on after the agents have already collapsed to the mean 

1 I -uf , opinion 0.5, one observes that, after a transient, the clus- 



12iV ■ iV V T ' * 



The asymptotic stationary solution correspond to 
Vt+T = Wt, a condition that immediately leads to Eq. 



ter spreads and the measured value of v agrees with the 
theoretical prediction ^. 

^.^ Aiming at further characterizing the system dynam- 
when plugged into pO|) . The above analysis also suggests 

7- TT- jS i~" u- u 4- 11 11 4- ics, we also studied the ease where, initially, agents share 

linear relation (I9|l , a choice which eventually allows us to carry ' j j ; & 

out the analytical calculation resulting in expression © (see ^he same belief Oq. The initial distribution of opinions is 

also the discussion in the appendixlA]) . Formally, eq. (0) applies therefore a Dirac delta f°{0) = S{0 - Oo)- Such condi- 

only for t = 0, when agents are populating the interval [0, 1] tion is a stationary solution, for any given Oo when the 

with a uniform distribution. During the subsequent evolution, death/birth process is inactivated. Conversely, when the 

the convergence still gives rise to a macroscopic triangular pat- death/birth applies, the system evolves toward a state, 

tern but, now, the associated triangle height Tc gets slightly characterized by a single cluster (localized, if p is small, 

reduced. At time t agents are still confined in the relevant in- uniformly spread over the allowed region as p > 1, see pre- 

terval [0, 1] and experience a certain degree of spreading, effect ^^^-^^ discussion), centered in O = 0.5 and with standard 

of the perturbation externally imposed. However, and espe- , . , . . i ^ m t^ • i i i ^i ^ ^i i- 

•^ I i- deviation given by Eq.[6J It is also observed that tire time 



cially for intermediate values of p, the progressive bunching 
opposes the birth/death disturb (which would tend to restore 
the i = variance) and drives an instantaneous reduction of 
Tc, as t grows. In the following, and to account for this self- 
consistent effect not captured by analytical framework, Tc is 
hence regarded as an effective parameter to be numerically 



needed by the system to complete the transition Tconv de- 
pends on the value of p and the critical affinity ac, see 
Fig. [51 A simple theoretical argument enables us to quan- 
titatively explain these findings. The initial distribution of 
opinions is modified after the first death/birth event as: 



adjusted: as commented below however, the best fit values of j^i/^\ "^ , N — Al 

/ (O) ^^+ ^ d{0-Oo). (10) 

Tc correlate well with direct measurements of the convergence 

(aggregation) time of the unperturbed system, a finding which The first term refers to the freshly injected actors, while 

a posteriori confirms the plausibility of eq. dSj). the second stands for the remaining Delta-distributed in- 



10 
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dividuals. Hence, the mean opinion value reads: 



Oi 



fHo)Odo 



M N ~ M 



2N 



N 



-Oo. 



(11) 



We can suppose that between the occurrence of two 
consecutive perturbations (separated by T iterations), the 
group average opinion docs not significantly change. No- 
tice that the probability of interaction of a newborn agent 
with another belonging to the main group is in fact pro- 
portional to M/N. Moreover several consecutive encoun- 
ters of this type are necessary to induce a macroscopic 
change of the averaged opinion. Under this hypothesis the 
next death/birth event makes the average opinion change 



of (fT5|) is in fact clearly demonstrated in Fig. [5^. Let us 
define the convergence time Tconv as the number of iter- 
ations needed to bring the average opinion e close to its 
asymptotic value 1/2. Solving Eq. ([T5|) for n and recalling 
that Tconv = nT yield: 



Tconv =T\0g 



\Oo 






(16) 



The above estimate is in excellent agreement with the nu- 
merical results reported in Fig. [SJs. 



as: 



- M N-M- 
O2 = \ Oi 

2N N 



(12) 



After n death/birth iterations, the opinion mean value 

reads : 

- M N-AI- 

On = ^ + ^^0„_, . (13) 

From Eq. ((T3)) one easily gets that the asymptotic equilib- 
rium is reached for Ooc ~ 0.5, as found in our numerical 
experiments; in fact the following relation is straightfor- 
wardly obtained: 

On-(^^) Oo + ^— (^^), (14) 
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being Oo the initial common believe. By setting a 
and /3 = -^ , the solution of Eq. P^ reads: 



1 - a" 

On = 13- + a"0„ , 

1 — a 



N-M 

N 



(15) 



whose asymptotic solution is given by 0„ — > Ooa = 0.5. 

Expression (jl5|) reproduces quite well the dynamics of 
the cluster mean, as seen in the simulations. The adequacy 



Fig. 5. a) The opinion as a function of time for the case where 
all the agents are initialized with the same opinion Oo = 0.1. 
The solid (red) line represent the numerical data while the 
segmented (blue) line represents the theoretical prediction from 
Eq. (|15|) . b) The transition time Tconv as a function of the death 
density p. Symbols stand for the results of the simulations and 
the solid line represents the prediction (|16p . Here A*' = 100, 
ac = AOc = 0.5, e = 0.001 and a = 0.07. 
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4 Conclusions to clarify this finding. As a general comment, we should 

emphasize that the effect of opening up the system to 
In this paper we have discussed the process of opinion 

external influences changes dramatically its intrinsic dy- 
making in an open group of interacting subjects. The 

namics revealing peculiar, potentially interesting, features 
model postulates the coupled dynamical evolution of both 

which deserves to be further explored, 
individuals' opinion and mutual affinity, according to the 

rules formulated in [7]. At variance with respect to the 

1 , r^ ^1 . . . ^ ^1 A Appendix 

toy — model [7j, the system is now open to contact with 

an external reservoir of potentially interacting candidates, ^j^j^ appendix is devoted to discussing a straightforward 

Every T iterations the M agents are instantaneously re- c^t^nsion of the above analysis to the case of the orig- 

placed by newborn actors, whose opinion and afhnity scores j^^j Oeffuant et al. model, which is made open via a 

are randomly generated according to a pre-assigned (here birth/death mechanism as outlined above. The interested 

uniform) probability distribution. The ratio p = M/T, ^^^^^^ ^^^ ^^^^^^^ ^ ^^^ ^ ^^^^.^ ^^^^^^^ ^^^^ ^^^ ^1^^^^ 

here termed departure density plays the role of a control ^^^^j specifications. We shall here solely recall that the 

parameter. The occurrence of a transition is found which D^ffuant's setting ^ = 1/2 and d = 1, is formally recov- 

separates between two macroscopically difl?erent regimes: ^^^^ ^y setting a, = into the aflinity model. 
For large values of the so-called social temperature the ^^ ^^^ ^j^^^^ Dcffuant's setting, assuming d = 1, the 

system collapses to a single cluster in opinion, while in the ^^^^^^^^ deviation of the opinion distribution decays as 

opposite regime a fragmented phase is detected. The role ^^ exponential function [H], namely: 

of p is elucidated and shown to enter in the critical thresh- 

v{t) = e-*/^'^w(0) , (17) 
old as a linear contribution. Two phenomena are then ad- 
dressed, with reference to the single clustered phase. On where Tc plays the role of a characteristic time. Dedicated 
the one side, the external perturbation, here being hypoth- numerical simulations, relative to the case study d = 1, 
esized to mimic a death/birth process, induces a spread- return the Tc = 191.52. 

ing of the final cluster. The associated variance is numeri- Assume now that every T-stcps the system opens: 

cally shown to depend on the density amount p, the tunc- Al agents are randomly removed. New actors enter the 

tional dependence being also analytically explained. On systems, their opinions being randomly sampled from a 

the other side, we also show that the birth/death events uniform distribution in the interval [0,1]. Let us denote 

imposed at a constant pace can produce the progressive p = M/T, the departure density. Furthermore, label with 

migration of a cluster, initially localized around a given v^°V^' the standard deviation of the opinions just before 

opinion value. A theoretical argument is also developed the insertion of the next pool of incoming agents. Hence, 
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in analogy with the preceding discussion, one can straight- 
forwardly write the recursive relation: 



^t+T 



M 
12N 



1-^] e-2^/ 



(18) 



whose asymptotic stationary solution corresponds to 



M 



12N 1 - (1 - M/N) 6-2^/^ 



(19) 



if T is small enough. This result can be compared to the 
result of direct numerical simulation, returning an excel- 
lent agreement, as displayed in Fig. [6] 



Remark 1 As a final remark, let us observe that in the 
general case, i.e. ac > 0, an equation formally analogous 
to (|18p can be derived, by invoking the correct exponen- 
tial ansatz (see main text). To obtain a closed analytical 
form for the asymptotic stationary standard deviation we 
however decided to resort to a linear approximation for 
the opinions convergence, as commented above. 

Clearly, when starting from a preformed cluster of opin- 
ions the injection of new actors determines an effective mi- 
gration of the mean, also ruled by Eq. (|13p in the original 
Deffuant et al. scheme. 



0,3 




Fig. 6. Final standard deviation, v, of the opinion cluster as 
a function of the departure density p. Each point results from 
averaging out 100 independent runs, each relative to A'' = 100 
agents. Symbols refer to numerical simulations performed with 
d = 1, M = 1/2 and different M : (o) M=l, (O) A/ = 2 and (□) 
Ad = 5. The solids line refers to the theoretical prediction (|19[l 
with Tc = 191.52, this latter being independently estimated as 
discussed in the main text. 
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